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ABSTRACT 

Diamond  has  been  grown  on  single  crystal  beryllium  oxide  by  hot-filament 
chemical  vapor  deposition.  The  diamond  particles  were  characterized  by 
Raman  spectroscopy,  and  the  morphology  and  microstructure  investigated 
by  scanning  electron  microscopy  (SEM)  and  plan-view  and  cross-sectional 
transmission  electron  microscopy  (TEM).  It  has  been  confirmed  that 
diamond  grew  epitaxially  on  the  basal  plane  of  BeO  with  the  epitaxial 

relationship  { 1 1  Udiamond/Zl 0001  IseO  and  <riO>diainond  rotated  by  less  dian6*’ 

with  respect  to  <ll20>BeO-  Diamond  was  also  grown  on  the  (1120)  prism 
plane  of  BeO.  A  selected  area  diffraction  pattern  obtained  from  plan-view 
specimens  indicated  the  presence  of  a  set  of  spots  that  corresponded  to  the 

[1010}  planes  of  hexagonal  diamond  or  to  the  ^(422)  forbidden  spots  of 

cubic  diamond.  During  diamond  deposition  on  (1 120)3eo  surfaces,  small 
particles  of  beryllium  carbide  (Be2C)  were  also  observed. 
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1.  INTRODUCTION 

Diamond  films  obtained  by  various  chemical  vapor  deposition  techniqnes 
are  almost  ail  polycrystalline  and  randomly  oriented.  However,  for 
electronic  applications  single  crystal  films  are  desirable.  Diamond  single 
crystal  films  with  large  areas  have  only  been  grown  snccessfuily  by 
homoepitaxy  on  single  crystal  diamond  substrates.  ^  Diamond 
heteroepitaxy  on  foreign  substrates,  e.g.,  Sr'^,  NF,  Cu^,  c-BN^-^®  and  0- 
has  been  limited  to  small  areas. 

Cubic  BN  has  a  similar  structure  to  diamond  and  close  lattice  match  (within 
2%).  However,  growth  of  large  single  crystal  c-BN  is  probably  as  difficult  as 
growing  large  diamond  single  crystals. 

Diamond  interfaces  with  c-BN  and  BeO  have  been  studied  theoretically. 
These  studies  show  that  BeO  is,  like  c-BN,  a  potentially  good  substrate  for 
heceroepitaxiai  growth  of  diamond.  The  Be-0  bond  length  in  BeO  is  0.165 
nm,  only  7%  greater  than  the  C-C  bond  length  of  0.154  nm  in  diamond. 
Also,  the  BeO/diamond  adhesion  energy  was  estimated  to  be  4.6  Jlm^, 
slightly  smaller  than  that  of  c-BN/diamond  (5.4 

One  advantage  of  BeO  is  that  high  quality  large  crystals  can  be  easily  grown 
by  a  flux  process. BeO  has  the  wunziie  crystal  structure^^  and  only  by 
heating  to  temperatures  above  2100*’C  does  it  transform  to  tetragonal  beta- 

beryllia.i3 

The  wurtzice  crystal  structure  is  based  on  a  h.c.p.  lattice  with  a  basis  of  two 
atoms  per  lattice  point;  Be  at  (1/3,  2/3,  0)  and  0  at  (1/3,  2/3,  3/8)  or  vice 
versa.  This  structure  (Table  1)  is  geometrically  similar  to  Lonsdaleiie 
(hexagonal  diamond)  and  thus  the  first  nearest  neighbor  positions  in  BeO 
are  similar  to  those  of  cubic  diamond.  Since  the  bonding  in  both  diamond 
and  BeO  is  dominated  by  nearest  neighbor  interactions,  BeO  is  a  particularly 
attractive  substrate  for  the  heceroepita.xial  growth  of  diamond. 

Very  different  situations  for  heteroepita.xy  are  encountered  on  the  different 

BeO  surfaces.  The  unreconstructed  prism  planes,  e.g.,  (1120)  and  (lOlO), 
present  a  situation  in  which  the  first  and  second  nearest  neighbors  are 
geometrically  similar  to  hexagonal  diamond.  These  surfaces  arc  comprised 
of  six-membered  rings  in  the  "boat"  conformation,  which  should 
presumably  favor  the  hetcroepitaxial  growth  of  hexagonal  diamond  (Fig.  1). 


The  unreconstructed  basal  plane  surface  of  BeO  is  comprised  of  six- 
membered  rings  in  the  "chair"  conformation.  The  first  nearest  neighbor 
interactions  are  ±erefore  the  same  for  a  layer  of  either  cubic  or  hexagonal 
diamond.  Also,  because  BeO  is  noa-centrosymmetric,  it  is  a  polar  crystal 
and  the  two  basal  planes  are  terminated  by  beryllium  and  oxygen.  We 
arbitrarily  choose  the  (0001)  surface  as  beryllium  terminated  and  the 

(0001)  as  oxygen  terminated.  In  an  atomic  hydrogen  environment,  the 
beryllium-terminated  basal  plane  is  expected  to  be  more  resistant  to  etching 
than  the  oxygen-terminated  plane.^^ 

n.  EXPERIMENTAL 

Hot-filament  assisted  chemical  vapor  deposidon  (HFCVD)  is  one  of  the  most 
popular  methods  to  grow  diamond  because  of  its  simplicity,  low-cost,  and 
ease  of  scaling  up  to  large  surfaces.  However,  metal  contamination  from 
the  filament  has  been  detected  within  the  deposited  film-*^  and  in  the 
interfacial  layers.-^---  Much  of  this  contamination  comes  from  the  inirial 
heating  up  of  the  filament  at  the  beginning  of  the  deposition,  perhaps  from 
volatile  metal  oxides  such  as  WO,.  To  reduce  contaminadon,  we  use  a 
movable  molybdenum  shutter  that  can  be  placed  between  the  filament  and 
the  substrate  during  stan-up.  Secondary  ion  mass  spectroscopy  (SIMS),  was 
used  to  detect  W  contamination  within  a  continuous  diamond  film  grown  on 
a  scratched  silicon  wafer.  No  W  contamination  was  detected  either  in  the 
bulk  of  the  film  or  at  the  interface  between  the  silicon  substrate  and  the 
film. 

Prior  to  start-up,  the  molybdenum  shutter  was  placed  in  between  the 
filaments  and  substrate.  A  flow  of  I(X)  seem  of  H2  at  20  Torr  was  staned 
and  the  temperature  of  four  filaments  connected  in  series  was  raised  to 
approximately  2000®C.  The  temperature  of  the  substrate  holder  during  this 
procedure  was  approximately  SOO^C.  Since  the  shutter  was  between  the 
filaments  and  the  substrate,  (he  substrate  surface  was  exposed  primarily  to 
molecular  hydrogen,  which  cleans  the  surface  of  nitrogen,  water, 
hydrocarbons,  or  any  other  atmospheric  contamination  adsorbed  by  the 
surface.  After  ten  minutes,  a  flow  of  0.5  seem  of  CH4  was  started.  The 
shutter  was  removed  5  minutes  later  and  the  distance  between  the  substrate 
holder  and  the  filaments  reduced  from  22  mm  [o5  mm.  The  sample 
temperature  was  estimated  by  a  thermocouple  at  one  side  of,  and  at  the 
same  level  as,  the  sample.  The  measured  temperature  was  850‘’C.  The 
substrate  holder  was  rotated  at  1/3  rpm  to  ensure  deposition  uniformity. 
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After  deposition,  the  samples  were  examined  by  conventional  TEM  in  a 
Phillips  CM20  scanning  transmission  eiecaon  microscope  (STEM)  operated 
at  200  kV.  High  resolution  transmission  electron  microscopy  (HRTE^0 
performed  in  a  JEOL  4(X)0EX  operated  at  400  kV  (point-to-point  resoludon 
of  -0.18  nm  at  Scherzer  defocus).  Micro  Raman  spectroscopy  was 
performed  with  a  Dilor  XY  modular  Raman  spectrometer,  which  has  a  spariai 
resolution  of  approximately  1  iim. 


m.  RESLXTS 

A.  Diamond  growth  on  (OOOl)BeO 

A  single  crystal  of  <0001>  oriented  BeO  was  provided  by  Brush  Wellman 
Inc.  After  verification  of  its  orientation  by  X-ray  Laue  reflection,  the  sample 
was  cut  into  small  pieces  (2x2x2  mm^).  The  (0001)  basal  plane  was 

selected  for  diamond  deposition  instead  of  the  (0001)  oxygen-terminated 
plane,  because  it  is  expected  to  be  more  resistant  to  attack  in  an  atomic 
hydrogen  environment.  We  determed  the  polarity  of  the  crystal  using  the 
pyroelectric  effect.^^  The  samples  were  polished  with  4000  grade  SiC  paper 
and  cleaned  ultrasonically  in  isopropanol.  SEM  and  Raman  spectroscopy 
were  performed  after  sample  removal  from  the  reactor.  Figure  2  shows  the 
hexagonal-shaped  crystals  obtained,  a  number  of  which  have  the  same 
orientation.  These  oriented  crystals  presented  a  strong  and  narrow  Raman 
line  at  1332  and  no  broad  peak  at  1500  cm*^  (Fig.  3). 

An  unreconstructed  step  along  a  <I0'10>  direction  on  the  basal  plane  of  a 

wurtzice  structure  defines  a  {2110}  prism  plane.  Assuming  that  diamond 
growth  on  a  stepped  (0001)  surface  of  BeO  lakes  place  by  nucleation  at  the 
step  followed  by  lateral  growth,  one  can  speculate  that  the  initial  deposit 
may  have  the  hexagonal  structure  (Fig.  4a).  An  unreconstructed  three  layer 
high  step  on  a  (OOOI)b40  su  rface  presents  the  same  geometrical 
configuration  as  boat-boat  bicyclodecane  (Fig.  4b)  proposed  as  a  model 
precursor  species  for  the  nucleation  of  diamond. Computer  simulation  of 
diamond  growth  from  nuclei  with  two  parallel  stacking  errors,  e.g.,  the 
boat-boat  bicyclodecane  shows  that  rapid  growth  of  a  thin  crystal  of 
Lonsdaleice  with  the  ABAB...  stacking  sequence  can  be  achieved  in  the  early 
stages  of  diamond  growth.-'^--^  After  the  elimination  of  the  surface  steps, 
subsequent  growth  on  this  initial  layer  would  presumably  be  cubic  diamond, 
which  is  more  stable  than  the  hexagonal  phase.  This  may  be  the  reason  that 
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the  oriented  crystals  chat  we  observed  in  this  study  are  hexagoaal  in  shape 
but  have  a  cubic  Raman  signal  at  1332  ctitL 

We  prepared  a  <1120>  cross-sectioaal  TEM  specimen  from  the  BcO  (0001) 
sample  for  conventional  and  high  resolution  TEM.  Figures  5a  and  5b 

respectively  show  a  cross-secrional  micrograph  and  the  corresponding  (llOl 
zone  axis  diffraction  pattern  of  cubic  diamond  with  the  electron  beam 

approximately  parallel  to  (1120lBeO-  To  reach  the  [llO]  zone  axis  of 

diamond  it  was  necessary  to  dlt  -d’  around  the  [0001]  direction  of  BeO. 

This  6’  misorientation  berween  the  (llO]  zone  axis  of  diamond  and  the 

[  1 1 20]  zone  axis  of  BeO  can  be  observed  in  Fig.  5c,  which  shows  the 

selected  area  diffraction  pattern  of  the  diamond/BeO  interface.  The  black 
arrows  in  this  figure  identify  the  III  diffraction  spots  of  diamond.  A  sec  of 
weak  spots  (white  arrows)  corresponding  to  0001  BeO  diffraction  spots 
planes  are  located  along  one  sec  of  III  diffraction  spots  of  diamond.  No 
other  structure  has  been  found  at  the  interface  between  the  diamond 
panicle  and  BeO. 

These  results  show  that  the  diamond  panicie  in  Fig.  5a  has  grown  epitaxially 
on  BeO  such  that  the  (ill)  planes  of  diamond  are  closely  parallel  to  the 

{0001}  planes  of  BeO  and  <ll0><iiamon<i  is  rotated  by  -6®  with  respect  to 

<ll20>BeO-  No  micros  true  tural  details  can  be  seen  in  Fig.  5a  because  the 
diamond  particle  is  too  thick. 

The  same  sample  was  further  treated  by  ion  milling  to  get  a  thinner 
interface  for  high  resolution  TEM.  Figure  6a  is  a  TEM  micrograph  of  the 
same  single  crystal  diamond  grain  shown  In  Fig.  5a.  At  the  top  left  comer  of 
this  micrograph,  a  micro-twin  band  can  be  observed  that  lies  on  (III) 
diamond  planes  parallel  to  the  (0001)  BeO  surface.  The  diffraction  pattern. 
Fig.  6b,  from  this  corner  clearly  shows  one  set  of  twin  spots  (white  arrows) 
located  at  one-third  of  the  distance  between  the  transmitted  beam  and  the 
primary  spots  from  the  unewinned  diamond;  this  is  typically  observed  in 
eVD  diamond  films. Figure  6c  shows  an  HRTEM  image  of  the  diamond. 

Since  the  diamond  grain  is  in  a  (1 10]  orientation,  two  sets  of  {III}  fringes 
with  0.206  nm  spacing  are  visible.  Note  that  micro-twins  are  only  present 
with  their  (111)  habit  planes  parallel  to  the  (0001)  surface  of  BcO.  The  6® 

misorientation  between  the  (ri0]diaraoBd  and  the  (U^OIbcO  zone  a.xis 
makes  it  impossible  to  obtain  a  good  HRTEM  image  of  BeO  and  diamond 
simultaneously,  and  thus  also  from  the  BeO/diamond  interface. 


Tilting  the  specimen  to  the  [11^]  zone  axis  of  BeO  (Fig.  6d)  made  i: 

possible  to  observe  the  (1100)  and  (0001)  fringes  of  BeO  with  0.234.  nm 

and  0.438  nm  spacings,  respectively.  The  (iH)  diamond  planes  visible  in 
Fig.  6d  form  an  angle  of  -70®  with  the  (0001)  planes  of  BeO.  This 
corresponds  to  the  angle  between  {Ill}  crystallographic  planes  of  cubic 
diamond  (70®32’)  and  shows  that  the  (111)  diamond  planes  are  parallel  to. 
the  BeO  (0001)  basal  plane. 

Similar  results  have  been  obtained  with  another  <1I20>  cross-sectional  TEM 
specimen  prepared  under  the  same  condirions.  Figure  7a  shows  a  cross- 
sectional  micrograph  of  a  diamond  particle  grown  on  the  basal  plane  of 
BeO.  In  this  micrograph  the  diamond  panicle  is  connected  to  the  BeO 
substrate  only  at  its  periphery  with  its  center  separated  from  ‘the  basal 
plane;  the  bright  contrast  at  the  interface  forms  an  angle  of  approximately 
5®.  Some  micro-twins  bands  parallel  to  the  BeO  surface  can  also  be 
observed  in  the  right  hand  side  of  the  diamond  panicle.  A  selected  area 
diffraction  pattern  (Fig.  7b)  was  obtained  from  the  right  hand  interface 
between  the  diamond  particle  and  the  BeO  substrate.  This  diffraction 
pattern  corresponds  to  a  <110>  zone  axis  of  cubic  diamond  with  the 

electron  beam  tilted  less  that  3®  with  respect  to  <1120>8eO-  A  set  of  spots 
(black  arrows),  corresponding  to  (0001  }b€0  planes,  forai  an  angle  of  -8® 
with  one  set  of  Ill  spots  of  cubic  diamond  (white  arrows).  We  believe  that 
at  the  early  stages  of  growth,  the  diamond  particle  has  grown  with  one  of 
its  (HI)  planes  parallel  to  the  (OOOl)BeO  surface.  During  the  process, 
internal  stress  produced  the  separation  between  the  diamond  particle  and 
the  surface  to  form  the  angle  of  8®. 

The  epitaxial  relationship  found  corresponds  to  the  ideal  superposition  of 
the  (111)  planes  of  diamond  on  top  of  the  (0001)  basal  planes  of  BeO  as 
show  in  the  model  presented  in  Fig.  8.  Here  it  is  assumed  that  the  carbon 
atoms  are  bonded  to  the  Be  atoms  at  the  interface. 


B.  Diamond  growth  on  (1120)B  e  0 

A  single  crystal  sample  of  BeO  with  a  (1120)  surface,  also  obtained  from 
Brush  Wellman  Inc.,  was  prepared  for  plan-view  TEM  prior  to  diamond 

deposition.  Several  3x3  mm-  sections  were  cut  from  the  original  (1120) 
oriented  sample  after  verification  of  its  orientation  by  X-ray  Laue  reflection. 
The  substrates  were  then  mechanically  thinned  to  a  thickness  of  100  m 


using  a  final  polish  with  4000  grade  SiC  paper.  They  were  then  dimpled  to  a 
thickness  of  30  .am  at  the  central  regioiL  Fmally,  ioo-milling  from  both 
sides  of  the  substrate  with  5  kV  Ar**  produced  a  hole  at  the  center  of  the 
sample.  At  this  stage  the  samples  were  introduced  to  the  CVD  reactor.  The 
deposition  conditions  were  the  same  as  those  used  for  (OOOl)BeO  substrate; 
only  the  deposition  dme  was  shoner  (4  hours)  in  order  to  study  the  inirial 
stages  of  nucieation.  Figure  9  shows  a  SEM  micrograph  of  one  of  the 
samples  prepared  in  this  way  after  Au  spunering  to  avoid  charging  of  the 
insulating  BeO  substrate.  Only  diamond  particles  grown  in  the  thinner 
regions  of  the  substrate  were  useful  for  TEM  studies.  This  micrograph  also 
shows  that  the  surface  of  the  sample  is  rough  and  stepped,  and  thus  it  is 
possible  that  some  of  the  diamond  parricles  did  not  grow  directly  on  the 

( 1  1  lO)  plane  of  BeO. 

A  dark-field  image  of  a  diamond  panicle  is  observed  in  Fig.  10a  with  the 

corresponding  selected  area  diffraction  (SAD)  pattern  (Fig.  10b)  obtained 
from  the  top  corner  of  the  diamond.  The  beam  direction  in  this  figure  is 

parallel  to  the  [1120]  zone  axis  of  BeO.  The  arrowed  spots  in  Fig.  10b 

identify  the  (1100),  (llOl),  (0001)  and  (0002)  planes  of  BeO  and  the  white 
ones  correspond  to  one  set  of  {111}  planes  of  cubic  diamond.  The  extra 

spots  observed  in  this  figure  come  from  double  diffraction  when  the 

original  electron  beam  is  diffracted  by  some  of  the  {111}  planes  of  the 
diamond  particle  on  top  of  BeO.  The  geometry  is  such  that  the  diffracted 
beams  from  the  diamond  particle  are  incident  on  some  of  the  BeO  planes  at 
the  Bragg  angle.  Effectively,  these  extra  spots  disappeared  when  the 
specimen  was  tilted  a  few  degrees  about  the  axis  connecting  the  111  spots 
of  cubic  diamond.  Under  these  conditions,  the  extra  spots  in  this 
diffraction  pattern  formed  a  similar  dark-field  image  as  the  one  shown  in 
Fig.  10a. 

The  same  diamond  particle  was  tilted  so  that  the  beam  direction  was 
parallel  to  the  [III]  zone  axis.  The  corresponding  diffraction  pattern  is 
shown  in  Fig.  lOc.  The  si.x  strong  spots  around  the  transmitted  spot 
correspond  to  {220}  planes  of  cubic  diamond  (d  =  0.126  nm).  Figure  lOd 
shows  a  dark-field  image  of  the  diamond  obtained  from  one  of  these  strong 

spots.  The  three  black-arrowed  spots  in  Fig.  10c  correspond  to  the  (I'lOO), 

(IlOO)  and  (0113)  planes  of  BeO.  In  addition,  six  weak  spots  (white 
arrows)  appear  around  the  direct  spot.  From  their  spacing  (d  =  0.219  nm), 

these  could  be  assigned  to  lOlO  spots  from  hexagonal  diamond  or  to  the 


^4.22)  forbidden  spots  from  cubic  diamond.  The  latter  could  occur 
because  of  incomplete  unit  cells  at  the  {III}  sanaces  of  diamond.-'^  As 

mentioned  before,  the  unreconstructed  (11^0)  prism  plane  presents  a 
"boat"  scrucnire  favorable  for  the  heteroepitaxiai  growth  of  hexagonal 
diamond.  In  this  case  a  buffer  layer  of  hexagonal  diamond  may  grow  at  the 
surface  of  BeO  before  a  cubic  diamond  structure  starts  growing.  Future 
work,  e.g.,  cross-secdonal  specimens,  will  be  necessary  to  confirm  the 
presence  of  this  buffer  layer. 

Small  panicles  of  beryllium  carbide  (BeiQ  ^cre  also  present  in  this  sample. 
Figure  11a  shows  a  SAD  pattern  obtained  from  a  BciC  particle.  The  020, 
220  and  200  diffraction  spots  of  BeiC  arc  indicated  by  black  arrows,  while 

the  white  arrows  identify  the  OOM,  IlOO  and  1101  diffraction  spots  of  BeO. 
This  indicates  that  the  electron  beam  is  parallel  to  the  [001]  direction  of 

Be-)C  and  the  [li20]  direction  of  BeO.  Figure  11b  shows  a  dark  field  image 
of  the  Be2C  panicle.  Note  that  some  small  BejC  panicles  also  light  up.  Tne 
structure  of  the  small  panicles  has  been  also  confirmed  by  the  HRTEM 

micrograph  shown  in  Fig.  11c  where  ^he  (lOlO)  and  (0001)  fringes  of  BeO 
with  0.234  nm  and  0.438  nm  spacing,  respectively,  are  observed.  The  BciC 

particle  shown  is  in  a  (110]  orientadon  and  two  secs  of  {1*11}  fringes  with 
0.247  n  m  spacing  are  visible.  Beryllium  carbide  has  the  andfluorite 
stnicture  (Table  1)  and  is  easily  hydrolyzed  by  contact  with  moist  air  to 
form  Be(OH)2  at  room  temperature.-*--’  Its  presence  indicates  that  BcjC 
can  be  stable  in  small  quantities,  or  when  it  is  protected  by  a  very  thin  layer 
of  another  material.  In  any  case,  its  formation  should  be  avoided  for 
diamond  heteroepitaxy  because  its  presence  is  not  favorable  for  this 
purpose . 


IV.  CONCLUSIONS 

The  conclusions  can  be  summarized  as  follows; 

(i)  Diamond  particles  were  grown  heteroepitaxially  on  the  basal  plane  of 
BeO  with  an  epitaxial  relationship  { 1 1 1 Idiamond //  (0001]b«o  ^nd 

<l*l0>diamond  rotated  by  less  than  6°  with  respect  to  <Il20>BeO-  Since  large 
crystals  of  BeO  can  be  grown,  this  material  is  a  potential  substrate  for 
growing  diamond  single  crystals,  or  highly  oriented  diamond  films. 


(ii)  HRTEM  image  and  SAD  panems  of  the  diamood/BeO  intenace  show 

chac,  under  the  condidons  adllzed,  ao  interfacial  phases  form  during  the 

deoosidon  of  diamond  on  BeO. 

* 

(iii)  No  epitaxial  reladonship  has  been  found  when  diamond  is  grown  on 
(H20)BeO- 

(iv)  Small  pardcies  of  BeiC  have  been  found  when  diamond  is  grown  on  the 
( 1  1 50)  prism  plane  of  BeO. 

(v)  SAD  pattern  of  diamond  grown  on  the  (1120)  prism  plane  of  BeO  shows 
the  presence  of  weak  spots  that  can  be  assigned  to  either  the  forbidden 

^422)  spots  of  cubic  diamond,  or  the  lOlO  diffraction  spots  of  hexagonal 
diamond. 
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HGURE  CAPTIONS 


Figure  1.  Cryscallographic  model  of  the  ideal  beteroepitaxial  growth  of 

hexagonal  diamond  on  the  (11^)  surface  of  BeO.  Grey,  white 
and  black  bails  are  respecnvely  carbon,  oxygen  and  beryllium 
atoms.  It  is  assumed  that  C  is  bonded  to  Be  at  the  interface. 

• 

Figure  2.  Hexagonal-shaped  diamonds,  some  showing  the  same  orientation, 
grown  on  the  (0001)  basal  plane  of  BeO. 

Figure  3.  Raman  signal  of  cubic  diamond  particles  shown  in  Fig.  2. 

Figure  4.  (a)  Crystallographic  model  representing  an  initial  deposit  of 
hexagonal  diamond  on  an  unreconstructed  three  layer  step  of 
(0001) BeO-  Grey,  white  and  black  balls  are  respecnvely  carbon, 
oxygen  and  beryllium  atoms. 

(b)  Boat-boat  bicyclodecane  with  two  parallel  stacking  errors. 
Positions  of  two  carbon  ad-atoms  are  shown  as  black  circles. 

Figure  5.  (a)  <11^0>bcO  cross-sccdonal  TEM  miCTOgraph  of  diamond  on 
BeO. 

(b)  (110]  zone  axis  diffracdon  pattern  of  the  diamond  single 
crystal  shown  in  Fig.  5a. 

(c)  [llOldiamond  SAD  pattern  from  the  BeO/diamond  interface 
shown  in  Fig.  5a. 

Figure  6.  (a)  TEM  micrograph  of  the  diamond  grain  shown  in  Fig.  5a 
after  further  ion  milling. 

(b)  Set  of  twin  spots  between  the  III  diamond  diffraction  spots 
corresponding  to  the  plane  parallel  to  the  (0001)  BeO  surface. 

(c)  HRTEM  image  of  the  [llO]  oriented  diamond,  {i'll}  fringes 
with  0.206  nm  spacing  are  visible  in  the  diamond. 


(d)  HRTEM  image  of  the  [11201  oriented  BeO  where  traces  of 


( 1  lOO)  and  (0001)  planes  with  0J234  /i/n  and  0.438  nm  spacing, 

respectively,  are  visible.  The  (ill)  diamond  planes  form  an  angle 
of  -70®  with  the  (0001)  planes  of  BeO. 

Figure  7.  Crystallographic  model  of  the  ideal  superposition  of  (111) 
diamond  on  top  of  the  (0001)  basal  planes  of  BeO.  Grey, 
white  and  blade  balls  are  respectively  carbon,  oxygen  and 
beryllium  atoms. 

Figure  8.  (a)  <1120>bcO  CTOSS-secdonal  TEM  micrograph  of  diamond  on 

BeO. 

(b)  [  1  lOldiamond  SAD  pattern  from  the  BeO/diamond  interface. 
Black  arrows  denote  BeO  000 1  difft-acdon  spots.  White  arrows 
denote  111  diffracrion  spots  from  cubic  diamond.  The  angle 
between  (0001)b«o  and  (lll)diamoad  is  approximately  8®. 


Figure  9.  SEM  micrograph  of  a  (1120)b«o  plan-view  TEM  sample  after 
diamond  deposidon. 

Figure  10.  (a)  Dark-field  image  of  diamond  grown  on  (1120)  surface  of 
BeO. 

(b)  SAD  pattern  obtained  with  the  electron  beam  parallel  to 

[1120]  zone  axis  of  BeO.  Black-arrowed  spots  identify  the 

( I  lOO),  (llOI),  (0001)  and  (0002)  planes  of  BeO;  the  white 
arrows  denote  correspond  to  one  set  of  {111}  planes  of  cubic 
diamond.  The  extra  spots  observed  in  this  figure  come  from 
double  diffraction. 

(c)  Diffraction  pattern  of  diamond  shown  in  Fig.  lOa  tilted  to  its 

[III]  zone  axis.  The  six  strong  spots  correspond  to  the  {220} 
planes  of  cubic  diamond.  The  three  black-arrowed  spots  are 
from  BeO  and  the  six  weak  spots  (white-arrowed)  can  be 

assigned  as  1010  of  hexagonal  diamond  or  the  ^422)  forbidden 
spots  of  cubic  diamond. 

(d)  Dark-field  image  of  the  diamond  obtained  with  one  of  the 
{220}  planes  shown  in  Fig.  lOc. 


Figure  II.  (a)  SAD  paaen  obtained  from  die  large  Be2C  panicie  siiowii  in 
Fig.  lib.  Blade  arrows  denote  die  020,  220  and  200  di&acdon 

spots  of  Bc2C  White  arrows  denote  the  0001,  IlOO  and  1 101 
diffraction  spots  of  BcO. 

(b)  Dark  field  image  of  BeoC  fonned  from  one  of  the  black- 
arrowed  spot  shown  in  Fig.  I  la. 

(c)  HRTEM  micrograph  of  a  small  Bc2C  particle  on  BeO  showing 

the  (1010)  and  (0001)  fringes  of  BeO  with  0J34.  nm  and  0.438 

nm  spacing  respecdvciy  and  two  sets  of  (ill)  fringes  of  Be2C 
with  0.247  nm  spacing. 


Figure  1.  Crystallographic  model  of  the  ideal  heteroepitaxial  growth  of 

hexagonal  diamond  on  the  (1 120)  surface  of  BeO.  Grey,  white 
and  black  balls  are  respectively  carbon,  oxygen  and  beryllium  atoms. 
It  is  assumed  here  that  C  is  bonded  to  Be  at  the  interface. 


Figure  2.  Hexagonal-shaped  diamonds,  some  showing  the  same  orientation, 
grown  on  the  (0001)  basal  plane  of  BeO. 


Figure  3. 


a) 


Figure  4. 


Raman  signal  of  cubic  diamond  particles  shown  in  Fig.  2. 


(a)  Crystallographic  model  representing  an  inirial  deposit  of 
hexagonal  diamond  on  an  unreconstructed  three  layer  step  of 
(OOOI)bcO.  Grey,  white  and  black  balls  are  respectively  carbon. 


oxygen  and  beryllium  atoms, 
(b)  Boat-boat  bicyclodecane 


Figure  5.  (a)  <1 1  20>bcO  cross-sectional  TEM  micrograph  of  diamond  on 
BcO. 


Figure  5.  (b)  [iTO]  zone  axis  diffraction  pattern  of  the  diamond  single 
crystal  shown  in  Fig.  5a. 


Figure  5.  (c)  [llOJdiamond  SAD  pattern  from  the  BeO/diamond  interface  shown 
in  Fig.  5a. 


Figure  6.  (a)  TEM  micrograph  of  the  diamond  grain  shown  in  Fig.  5a  after 
further  ion  milling. 


Figure  6.  (b)  Set  of  twin  spots  between  the  1 1 1  diamond  diffraction  spots 
corresponding  to  the  plane  parallel  to  the  (0001)  BeO  surface. 


uv 


f-ijure  6.  (d  )  HRTEM  image  of  :he  [1  l!0i  onented  BeO  where  traces  of 

( !  [(;())  and  (ChOl)  planes  winh  0.234  nm  and  0.43S  nrn  spacing, 

rv.j^pev.t:ve!y,  are  visible.  The  (111)  diamond  planes  form  an  an^Ie  of 
-.'O'"  With  the  (0001)  planes  of  BeO. 


Figure  7.  Crystallographic  model  of  the  ideal  superposition  of  (1 1 1)  diamond 
on  top  of  the  (0001)  basal  planes  of  BeO.  Grey,  white  and  black 
balls  are  respectively  carbon,  oxygen  and  beryllium  atoms. 


Figure  8.  (a)  <ll20>BeO  cross-sectional  TEM  micrograph  of  diamond  on  BeO. 


Figure  8.  (b)  [llO]diamond  SAD  pattern  from  the  BeO/diamond  interface.  Black 
arrows  denote  BeO  0001  diffraction  spots.  White  arrows  denote  111 
diffraction  spots  from  cubic  diamond.  The  angle  between  (OOOl)BeO 
and  (lll)diamor\d  is  approximatly  8-. 


Figure  10.  (a)  Dark-field  image  of  diamond  grown  on  (1 120)  surface  of  BeO. 


Figure  10.  (b)  SAD  pattern  obtained  with  the  electron  beam  parallel  to  [1120] 

zone  axis  of  BcO.  Black-arrowed  spots  identify  the  (1 100),  (1  lOl), 
(0001)  and  (0002)  planes  of  BeO;  the  white  ones  correspond  to  one 
set  of  ( 1 1 1 )  planes  of  cubic  diamond.  The  extra  spots  observed  in 
this  figure  come  from  double  diffraction. 


Figure  10.  (c)  Diffraction  pattern  of  diamond  shown  in  Fig.  10a  tilted  to  its 

[111]  zone  axis.  The  six  strong  spots  correspond  to  the  {220}  planes 
of  cubic  diamond.  The  three  black-arrowed  spots  are  from  BeO  and 

the  six  weak  spots  (white-arrowed)  can  be  assigned  as  lOlO  of 
hexagonal  diamond  or  |(422)  forbidden  spots  of  cubic  diamond. 


Figure  10.  (d)  Dark-field  image  of  the  diamond  obtained  with  one  of  the  {2201 
planes  shown  in  Fig.  10c. 


Figure  11.  (a)  SAD  pattern  obtained  from  the  large  Be2C  particle  shown  in 
Fig.l  lb.  Black  arrows  denote  the  020, 220  and  200  diffraction 

spots  of  Be2C.  White  arrows  denote  the  0001,  (IlOO)  and  (IlOl) 
diffraction  spots  of  BeO. 


Figure  11.  (b)  Dark  field  image  of  Be2C  formed  from  one  of  the  black-arrowed 
spot  shown  in  Fig.  1  la. 
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Figure  1 1.  (c)  HRTEM  micrograph  of  a  small  Be2C  particle  showing  the  (1010) 
and  (0001)  fringes  of  BeO  with  0.234  nm  and  0.438  nm  spacing 

respectively  and  two  sets  of  (111)  fringes  of  Be2C  with  0.247  nm 
spacing. 
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Electronic  5tructure  calculations  predict  that  hypothetical  ordered  mixed  crystals  of  c-BN  and 
dianiond-C  representative  of  the  expecteti  short-range  order  in  the  alloys  shotr  a  very  pronounced 
hand-gap  bowing.  This  anonialons  behavior  is  believed  to  be  related  to  the  staggeretl  hand  Gne-up  at 
the  correspontling  heterojnnction.  The  miscibility  phase  diagram  of  this  system  is  estimatetl  u,«ing  a 
sinipie  psendohinary  'regular  solution"  model  for  the  short-range  order  and  the  energies  of  fonitation 
of  the  representative  orderetl  compounds  calculated  from  first-principles.  It  indicates  that  there  is 
only  very  limited  mutual  solubility  in  the  solid  stale.  The  bandgaps  and  energies  of  formation  of 
the  disordereil  ternary  alloys  are  estimated  by  means  of  cluster  expansions. 

P.AfS:  71.  >.5.Tn.  61.5.5Hg 


I.  INTRODUCTION 

Diatiioiid  and  c-BN  have  recently  received  con.sider3ble  attention  becau-se  of  their  promising  nialeriaU  properties 
for  alira.sives.  lieat.sink.s.  protective  coatings  and  even  wide-band-gap  semiconductor  applications.  This  interest  stems 
from  tlie  extreme  values  of  these  two  materials'  properties  such  as  hardness,  thermal  conductivity,  elastic  coiistant.«. 
ami  l>aiul  gaps  (!].  The  recent  interest  has  surged  as  a  result  of  the  increased  availability  of  these  difficitit  to  obtain 
materials.  This  in  turn  resulted  from  the  development  of  a  number  of  novel  growth  techniques  among  which  are 
energetically  enlia'ticed  (plasma,  hot-filament,  microwave)  chemical  vapor  deposition  (CVD)  (2],  laser  ablation  (3]  ami 
ion  beam  techniques  (4).  A  basic  difficulty  common  to  both  materials  is  that  they  are  tliennodynaniically  only  staldi- 
at  Iiigli  pressures,  the  low-pressure  ground  states  being  the  layered  hexagonal  phases:  graphite  and  h-BN. 

Mixetl  (.’-BN  layered  hexagonal  phases  have  been  synthesized  by  a  number  of  groups  [5-7]  and  studietl  theoretically 
by  Liu  el  al.  [8j.  Batlziaii  [-5]  lia.s  .synthesize^I  mixed  crystals  of  diamond  and  r-BN  by  a  high-prrssiirr  high  lemptraliire 
plia«e  trnii.sformation  technique  starting  from  the  layered  he.xagonal  mixed  crystals.  Attempts  to  grow  riibic  mixed 
pha.ses  l\v  (,‘VD  Itave  so  far  failed. 

In  this  paper,  the  results  of  electronic  structure  calculations  of  the  cubic  (BN),C2(i-v»  sy-^feni  are  presented.  Tlii.' 
study  complemeiils  a  previous  .study  of  the  diamond/c-BN  lieierojunclion  [9j.  Our  calculations  address  the  quest  ion  of 
tin’  band  ga(>  behavior,  the  energy  of  formation  and  the  phase  stability.  Rather  restrictive  but  plau-sihle  assumptions 
an-  made  about  the  sliort-range  order  in  this  alloy  .system. 

Tile  paper  is  organized  as  follows.  Sec.  1!  describes  the  computational  approach.  .Sec.  Ill  desrr  bes  the  assumptions 
about  the  .sort-range  order.  In  section  IV,  we  first  dlscass  the  results  for  the  band  gaps  of  the  ordered  coiupouiuls, 
next  I  lie  plia.se  diagram  of  the  alloys,  and  finally,  the  properties  of  the  disordered  alloys.  The  conclusions  of  this  work 
are  summarized  in  .Sec.V. 


II.  (.  O.MPUT.^TION.^L  .METHOD 


Tli»’  results  presented  here  were  obtained  using  Andersen’s  linear  muffin-tin  orbital  method  (LMTO)  in  the  atomic 
sphere  a|ipru.\imatioii  ( A.SA )  itirliidiiig  the  .so-called  combined  correction  [lUj.  Tlie  Koliii-.Sliam  local  density  fiinrt  ioiial 
theory  (1 1]  is  used  willi  the  exchange-correlatioii  paramefrizafion  of  Hedin  and  Lumiqvist  (r,fj. 

The  local  density  approximation  (LDA)  is  well  known  to  underestimate  the  band  ga|is  becaase  the  Kohii-.Sham 
eigenvalues  are  not  strictly  speaking  quasipartirle  excitation  energies.  In  this  work,  we  estimate  the  qiia.siparlicle  .self- 
energy  correction  using  an  expression  derived  by  Bechstedt  and  Del  .Sole  [l.3j.  The  latter  Is  ba.sed  on  an  approximate 
Ireatiiieiit  of  Hediii's  GW  approximation  [14].  In  Becbstedf  and  Del  Sole’s  treatment  (l3j.  local  field  efTecIs  are 
neglected  aiid  matrix  elements  are  calculated  with  a  simple  tight-binding  approximation  to  the  LDA  wave  fuiictioii.s. 
Ill  the  final  closed  expres.sion  of  the  correction. 
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only  wfll-kiio'vii  (juaiitities  app«far  such  as  the  dielectric  constant  the  Thomas- Feniii  screening  wave  vector  'ur- 
ami  1%//  =  [( 1  —  np)r^/2  -(•  ( 1  +  up)rg/2]  with  r^  and  rg  effecive  radii  of  cation  and  anion  Slater  orbitals  given 
hy  ’■-(/£»  =  (l/■l”(1.7  T  0.05|Z^  -  Za\)  with  a  the  cubic  lattice  constant  and  the  bond  polariaahilities  op  has«t  on 
Harrison's  universal  tight-binding  theory  [15]. 

Our  approach  to  studying  the  disordered  alloys  involves  two  parts.  The  first  is  the  calculation  of  properties  of 
ordered  compounds  representative  of  the  assumed  short-range  order.  In  the  second  part,  the  statistically  average<l 
properties  of  ihe  disordered  alloys  are  determined  by  a  cluster  e.^pansion  whose  coefficients  are  derived  by  a  mapt'ing 
of  the  cluster  expansion  to  Ihe  ordered  crystal  results.  This  is  a  generalization  of  the  Connoly- Williams  approach] Hi] . 
In  I  he  follow  ing  section,  we  discu.ss  our  assumptions  about  the  short-range  order  and  the  choice  of  ordered  stnictures. 

III.  MODEL  OF  SHORT-RA.MGE  ORDER 

The  short-range  order  in  leterahedrally  bonded  heterovaleut  ternary  alloys,  of  which  the  presently  studieil  (  '-BN 
.system  is  a  particular  ca.se.  has  been  discu.ssed  vigoroasiy  in  a  .set  of  papers  mostly  in  the  mifidle  ttp’s.  E.xlettsiv>* 
experimental  and  theoretical  studies  were  devoted  to  (GaAslj-Ge^u.,,.  (Ga.Sb)j.Ge-ju_,,  [17-215].  A  recent  overview 
of  the  fiehl  can  be  found  in  Ref.  [2(5]. 

In  tlie  dilute  limit  of  the  111-V  component  in  the  alloy,  one  may  expect  isolated  group-III  and  grotip-V  impurities  in 
a  diamond  lattice.  That  is,  one  can  expect  group-111  and  grouivV  elements  to  occur  on  either  of  the  two  fcc-sublattices 
of  the  diamoiul  lattice  with  ecpial  probability.  At  a  higher  concentration  of  the  Ill-V  component  one  may  expect  a 
transition  to  a  nominally  zinc-blende  lattice  in  which  the  group-111  element  has  a  statistical  preference  for  one  of 
the  sublatlices  and  the  group- \‘  element  for  the  other.  The  simplest  way  to  view  this  transition  Ls  as  a  percolation 
problem  [H?].  More  sophisticated  treatments  will  include  energetic  and/or  .surface  growth  considerations. 

From  an  energetic  point  of  view,  one  may  expect  that  the  group-III  and  group-V  elemenLs  would  tend  to  be  nearest 
neiglibors.  This  is  becau.se  the  III-IV  and  IV-V  (AZ  =  il)  bonds  are  under-  and  oversaturated  respectively  by 
1/4  electron/bond.  One  thus  can  gain  up  to  —1/4  of  the  bandgap  fy  by  charge  transfer  from  one  to  the  other. 
This  compensation  leads  to  charged  ‘ijnor  and  acctpior  bonds  whicli  costs  electrostatic  energy.  As  a  result,  the 
compensation  is  mast  effective  for  nearest  neighbor  III-IV  and  IV’-V  bonds  in  which  case  the  electrostatic  energy  cost 
Is  minimized.  These  compensation  effects  between  over-and  uudersaturated  bonds  have  been  discassed  by  Dandrea 
et  al.  [27]  and  Lambrecht  et  al.  (28). 

For  .similar  electrastatic  reasons.  AZ  =  ±2  bouds  ( V-V  and  IlI-lII)  are  unlikely  to  occur.  These  basically  correspond 
to  antisite  defects,  which  have  a  large  energy  of  formation  (21)].  The  same  conclusion  is  reached  from  a  study  of  inversion 
domain  boundaries  in  GaAs  [28].  In  the  ca.se  of  the  Ge-based  alloy.s.  there  has  been  a  controversy  over  the  pres«'nre 
or  ab.s*»iice  of  the  AZ  =  ±2  so-called  '‘wrong”  bond.s  [20.22].  It  was  shown  by  Holloway  and  Davis  [2‘2j  that  their 
presence  in  any  significant  amount  would  lead  to  a  closure  of  the  gap  becaii.se  of  the  associated  deep  levels.  This 
wouhl  lie  in  disagreement  with  the  e.xperimental  data  on  the  bandgap  beliavior  in  these  alloys  and  this  was  ii-serl  as 
an  argument  for  their  alisence. 

Osorio  et  al.  [2G]  included  the  above  mentioned  electrostatic  effects  explicitly  in  a  statistical  ihermcKlyiiauiirs 
trcntmi’iil  of  tli^s*’  alloy.s.  They  concluded  that  the  electrostatic  effects  lend  to  stabilize  the  zinc-blende  ordering  over 
the  diamoiul  structure.  Belo'v  tlie  melting  temperature,  however,  botli  pha.ses  were  found  to  he  unstable  towards 
phase  separation  over  most  of  the  composition  range.  The  energetic  arguments  against  the  AZ  =  ±2  bonds  and  for 
tlie  “olnstering"  of  AZ  =  ±I  bonds  are  even  stronger  for  compounds  with  larger  ionicity  and  thus  are  more  compelling 
for  c-BN  Ihnn  for  the  GaAs  and  GaSh  cases  studied  earlier. 

From  Bnilzi.iu's  .\-ray  diffraction  studies  [5].  it  appears  that  at  lea.'st  partial  zinc-blende  ordering  can  be  ascimied 
for  the  solid  solutions  of  ctduc  C-BN  prepared  by  him.  That  is.  the  B  nud  N  ions  maintain  a  preferential  ocrnpaiion 
on  spp.Trale  frc  snldaltices  of  wiiirli  a  certain  percentage  are  replaced  by  carbon  atoms  in  Ihe  mixed  rryslal. 

In  view  of  the  above,  we  a.ssume  that  the  short-range  ordering  is  characterized  hy  (1)  absence  of  B-B  and  .N- 
N  boinls.  ami  (2)  nearest  neighbor  positioning  of  B  and  N  leading  to  local  compen.salion  of  B-C  ami  N-C  bonds. 
These  plausible  assutiiptioiis  correspond  ba.sically  to  a.s$itmiug  an  altiiosi  perfect  local  chtirgf  nmlralilii  ami  lof^al 
.^Inirliioiiiftn/.  Obvioitsly.  this  cannot  be  strictly  applicable  in  Ihe  dilute  limit  of  .small  r  in  (BN),C;ii-i-).  We  choose 
here  not  to  treat  the  controversial  subject  of  Ihe  zinc-hleude-to-dianioiKl  transition.  Rallier.  we  focas  on  the  coars-T 
prolilem  of  the  miscibility  and  the  properties  of  the  most  plau-sible  .short-range  order  structure  over  most  of  the 
composition  range  (cxclndiiig  r  <  1).  which,  as  argued  above.  ha.s  the  zinc-bleiide  .structure. 

As  a  starling  point,  we  study  hy|>olhetical  orderetl  structures  ba.sed  on  tlie  zinr-blende  struclnre.  The  models  we 
investigated  correspond  for  r  =  0.5  to  the  Li,  and  for  r  =  0.25  and  r  =  0.75  to  the  ii,  frr-.snperlallires  of  BN 
ami  (  ■;  "molernles"  oriented  along  the  (111)  direction.  Of  course,  this  does  not  imply  that  we  altlarh  a  meaning 
to  these  molecules  as  independent  or  weakly  interacting  entities.  The  .sy.stem  obviously  is  not  a  molecular  rrysial. 
This  fhoife  corresponds  to  the  Coiinolly-Williaim  |C\V)  [16]  structures  with  »  £  {0.  l.'i.3.4}  on  one  of  tin- 
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fcc-siil>lafticP5  In  '•acli  of  thww.  then*  occurs  exactly  one  type  of  the  BnC^-n  fcc-nearest  neighbor  teiralieUra.  Thteu* 
leirahetlra  form  the  basis  of  the  cluster  expansions  used  for  the  statistical  tratnient  of  the  disordered  alloys.  For  the 
purpose  of  statistically  averaging  the  properties  of  the  alloys,  there  is  a  one-to-one  correspondence  betwe.Mi  the  five 
liasic  tetrahetlra  and  the  five  CW-structures. 


IV.  RESULTS 

A.  Band  gaps  of  ordered  coiiiponmLo 

We  start  with  a  ili.sciission  of  the  bai  '  .ps  of  the  pure  materials.  The  bancigap  of  c-BN  has  been  somewhat 
controversial  in  the  pa.st  becau.se  of  the  difficulty  of  obtaining  single  crystals  of  adequate  size.  A  value  of  6.4  eV'  as 
obtained  by  ab.corptioii  measurements  by  Chrenko  [-301  on  high-pressure  .synthesizetl  cyrstals  [31].  Recently.  .Mishima 
[3‘2j  mea.siired  the  optical  ah.sorptiou  edge  to  be  61  ±0.2  eV.  Our  previous  calriilatioiLs  yieideil  a  band  gap  of  1.4(5  eV 
within  LDA  and  (5.4  eV  iiirludiiig  the  Becltsfedf-Del  .Sole  correction.  The  band  gap  of  diamond  obiaiiierl  by  the  same 
computational  approach  was  -5.6  eV  (4.1  eV  in  LDA)  in  fair  agreement  with  the  experimental  value  of  5.5  eV.  From 
our  ami  their  experieiic**.  the  Bechstedt-Del  Sole  approach  is  usually  found  to  give  a  sliglit  overestimate  with  errors 
of  the  order  of  a  few  01  eV  at  most. 

Ba.sed  on  the  fart  that  the  band  gap  of  c-BN  is  larger  than  that  of  diamond,  one  might  e.vpect  that  mixfd  diamoini- 
c-BN  cry.stals  would  have  baud  gaps  larger  than  that  of  ciiamond.  The  calculated  baud  gaps  of  ordered  compounds 
are  shown  in  Fig.  1.  The  LDA  values  are  indicated  by  squares,  the  corrected  values  [33]  by  diamonds.  Contrary 
to  simple  expectations,  one  .sees  that  the  gap  initially  shows  a  dramatic  decrease  upon  alloying  c-BN  into  diamond. 
Even  for  the  75  %  c-BN  compound,  the  gap  is  smaller  than  that  of  diamond. 

We  believe  that  the  rea.son  for  this  anomalously  strong  bowing  of  the  bandgap  can  be  understood  in  terms  of  the 
alignm>‘nt  of  the  energy  bands  in  c-BN  and  diamond.  As  was  shown  in  our  earlier  work  on  the  heterojunctioiis  betwefii 
these  two  solitls.  the  baud  offset  is  of  the  so-called  type-II.  (i.e.  staggered  type),  as  shown  in  Fig.  2.  This  means 
that  the  valence-I>aud  maximum  of  the  composed  system  (muced  crystal  or  heterojuiictiou)  is  mainly  C-like  while  the 
coiidurtioii-baiid  minimum  is  mainly  BN-like  with,  in  fact,  a  significant  B-componenl  because  of  the  cation  nature  of 
the  latter.  This  implies  that  the  effective  band  gap  of  the  system  b  lower  than  that  of  either  of  tlie  two  pure  materials. 
Tile  effect,  in  fact,  is  slightly  more  pronounced  in  the  ii,  structure  studied  here,  which  is  a  1  -H  (001)  superialtire, 
than  it  is  in  larger  period  superlattices  approaching  the  semi-infinite  heterojunctioii  [0].  This  effect  clearly  may  be 
expectetl  to  occur  in  other  alloy  systems  in  which  the  heterojunctioiis  between  the  two  constituent  srinicomhirtors  is 
of  type  n.  e.g,  Sit  '/BP  [34]. 


B.  Phase  slabtlilr 

As  .ilrpaily  m*'ntioned.  it  is  not  clear  that  this  sy.steiii  is  a  true  solid  solution  in  the  tliermodyiiamir.seu.se.  First  of 
all.  tlii.s  appears  not  to  be  the  case  for  the  other  heterovaleni  ternary  alloy  systems,  and  secoiully  one  may  e.'cperl  pli.i.si- 
.sf'paralion  on  tin-  ba.sis  of  the  high  interface  energy  of  formatioii  that  we  previou.sly  calculated  for  the  heterojunction 
[!)].  We  discuss  the  miscibility  phase  diagram  here  within  a  very  simple  approximation. 

First,  ill  accordance  with  the  as.sumptions  about  short-range  order  discassed  in  Sec.  III.  we  treat  the  alloy-problem 
as  a  pseuilobinary  system  of  BN  and  C’  ■molecules".  Secondly,  the  total  energy  of  the  solid  solution  is  expanded  in 
.1  f  luster  expansion  using  the  above  described  fcr-tetrahedroii  appro.ximation.  The  cluster  energies  are  identified  witli 
llie  eiu  rgies  of  format  inn  of  the  corresponding  C '^^’-sfr^lctures  as  explainer!  in  Sec.  III. 

For  sitnplirity,  we  ignore  the  small  lattice  mismatch  (~  1  %)  between  diamond  and  c-B.N.  Thus  we  only  need  tin- 
energy  of  formation  at  their  equilil'riiim  volume.  These  are  found  from  our  first-principles  total  energy  ralrulalions 
to  follow  a  paral>ola  very  closely,  that  Is 

A//,.(r)  =  4Af/-..r(l -r)  (2) 

with  AWi  the  energy  of  the  Li„  compound  equal  to  U.45e\7atom.  The  subscript  “o“  staiid.s  for  “ordered"  com|>outid 
as  opposerl  to  the  disordered  alloys.  While  in  general  the  CW-cluster  expansion  is  equivalent  to  a  ffntralizti  Isitig 
spill  model  iiicliiiling  triplet  and  tetrahedron  interactions,  the  special  case  of  a  parabolic  dependence  on  r  Is  well 
known  to  correspond  to  the  Isiiig  model  with  only  nearest  neighbor  pair  interactions. 

H  =  (3) 

(’/) 

III  tlie  present  case.  7;  >  0  since  like  “molecules"  prefer  each  other  over  unlike  ones  as  nearest  neighbors.  This 
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is  imlicat*^!  by  the  positive  energy  of  formation  of  the  compounds.  This  means  that  the  alloy  statistics  in  this 
approximation  is  mapped  on  that  of  the  well  known  ferromagnetic  Ising  model. 

\V(*  liriefly  recall  the  *regidar  solution"  or  meaii-fiehl  treatment  of  this  problem.  In  this  model,  the  entropy  of 
mixing  is  approximated  by  the  point-appro.ximation  in  which  S  =  -t(rlnr  -  (1  -  /)  ln(l  -  r)]  with  k  the  Boitzman 
constant.  The  free  energy  becomes  F{t.  T)  =  A//,(r)  -  TS{r).  At  low  temperature,  this  fuiirtion  shows  a  maximum 
at  r  =  tj.')  and  two  minima  symmetrically  located  on  either  side  as  shown  in  Fig.  -I.  With  increasing  temperature,  the 
two  ininitna  converge  towards  the  inaximnm  until  above  a  certain  critical  temperature.  they  merge  together  in 

a  single  Miininuiin  at  >0  %.  Tlie  locus  of  these  minima  in  the  r  -T  plane  describes  the  mi.sfibility  gap  (hinodal  line) 
and  i.s  given  by 

k-TI^Hi  =  (8r-.4)/(lnr-ln(l-r)].  W 

The  region  where  d'  F/dr-  <  U  is  unstable  and  is  boundeil  by  the  spinodai  line 

i-r/A//;  =  8r(l -j).  l-^)) 

The  critical  temperature  T,\ig  —  1-^Hilk  as  7700  K. 

Tlie  ferromagnetic  Ising  model  has  been  treated  by  means  of  more  sophisticated  statistical  methods,  such  as 
Kiknrlii  s  clii.ster  variation  method  (CVM)  (.3‘j..36].  With  increasing  size  of  the  clusters  taken  into  arcotint  in  the  en¬ 
tropy  expression,  one  olitains  a  decreasing  correction  factor  of  O-SS-S-O.^SO  for  Tma  relative  to  that  for  the  regular  solu¬ 
tion  approximation.  The  most  accurate  value  known  for  the  critical  temperature  was  obtained  by  a  high-temperature 
expansion  [-37]  wliicli  gives  the  critical  temperature  as  0.81627  times  tlie  critical  temperature  of  the  regular  .sobitinn 
model.  Even  with  this  correction  factor.  Tsto  *=  6300  K. 

This  is  much  higher  than  the  melting  temperature,  which  was  estimated  to  be  —3500  K  for  c-BN  and  —4000  K  for 
diaiiioiul  at  —10  GPa(l].  A  tentative  phase  diagram  including  the  renormalization  of  the  critical  temperature  and  the 
position  of  the  liquidus  line  is  shown  in  Fig.  4(a).  As  can  be  seen,  our  model  predirts  very  limited  mutual  sotiibility  of 
c-BN  and  diamond  in  the  solid  state,  in  apparent  conSict  with  Badzians  claim  of  having  synthesized  solid  solutions 
of  compo-sitions  0.16  <  ^bn  <  0.6. 

We  next,  consider  possible  improvements  to  the  simplified  trealnieut  given  above.  Including  the  lattice  misniatch  of 
about  1  %  requires  expanding  diamond  and  compressing  c-BN  to  the  common  lattice  constant  of  the  murti  crystal 
which  will  cost  an  additional  energy  A£r-  This  should  be  followed  by  a  rela.xation  of  the  bond  lengths  in  Ih**  alloy, 
releasing  an  energy  A£r.  These  two  effects  partially  cancel  each  other.  We  estimate  that  this  may  change  the  energies 
of  formation  and  Mins  the  temperature  scale  by  at  most  a  few  percent.  This  would  thus  not  alter  our  ronrhision  lliat 
the  inisribilil.y  critical  temperature  lies  well  above  the  melting  temperature. 

A  more  important  consideration  might  be  the  inclusion  of  long-range  correlations.  In  some  coiiinioii-aiiion  psHU- 
ilobiiiary  .systems,  this  lia.s  been  .shown  to  lead  to  lowering  of  the  mi.scibility  rritiral  temperature  by  as  ninch  a.s  60 
%  [38].  A  reiiornialization  of  this  magnitude  would  bring  T.ug  down  to  the  order  of  the  melting  temperature.  In 
particular,  llie  fourth  nearest  neiglibor  pair  interaction,  which  is  between  cations  linked  along  the  bond-chains  in  tln- 
(110)  dirt'Clion.  was  found  to  be  quite  important  in  Ferreira  et  al.’s  [.38)  work.  .Since  this  has  been  foiiml  to  b**  a 
rat  her  gener.Tl  feature  [38.  ^O).  a  .similar  effect  may  well  occur  in  the  present  .sy.stem.  For  the  purpose  of  iIln.slr.Tt  ing 
the  following  discu.s.sion.  we  show  in  Fig.  4(h)  another  tentative  plia.sediagraiii  which  includes  a  60  %  renofiiialization 
due  to  these  effects. 

A  complete  treatment  of  long-range  effect.s  is  beyomi  the  scope  of  the  present  paper.  Snrii  a  treatment  .should 
prob.ibly  be  ba.'.ed  on  extensions  of  the  Blnme-Eiiiery-Griffiflis  .5=1  spin  Hamiltonian  [40]  wliicli  was  first  applied 
t'j  Mie  problem  of  In  lerovaleiit  ternary  (III-\')-IV  allo.vs  by  Newman  el  al.  [20].  The  generic  plia.se  diagram  of  this 
model  flamilloniaii  is  characterized  by  the  presence  of  a  nii.scibilily  gap  ami  an  order-di.sorder  Iransiiion  b«-lween  tin- 
zinc-blende  and  diamond  pha.ses  in  the  .solid  solution  region  alwe  the  mi.scibility  gap  which  ends  at  the  binod.il  line 
in  a  tricritiral  point  We  have  thus  tentatively  indicated  a  ziiic-bleiide-to-diamond  transilioii  in  Fig.  4(b)  with  the 
understanding  tliat  tlie  position  of  tlie  transition  is.  of  course,  not  known  at  this  point.  a.s  is  indicated  by  the  (|neslioii 
mark.  We  tentatively  |)lace<l  Mie  line  close  to  the  diamond  limit  becau.se  Badzian  s  [5)  resnlls  indicate  that  al  lea.*:! 
partial  zinc-blemie  ordering  was  present  in  his  samples  down  to  c  =  0.1-').  We  empha-size  that  the  form  of  the  plm.se 
diagram  shown  here  is  very  scliem.atic.  .See  for  eg.  the  papers  of  Newman  et  al.  [2Uj  and  Osorio  et  al.  (26)  for  realistic 
plia.se  diagrams  of  tiie  ,5  =  1  Hamiltonian  for  various  values  of  the  relevant  paranielers. 

Ill  the  .situation  depicted  by  Fig.  4(b).  one  may  imagine  quenching  the  .sy.stem  from  above  the  miscibility  gap 
to  lower  leiiiperainres.  Above  the  miscibility  gap.  but  belo'v  the  melting  line,  the  system  would  either  have  the 
diamond  or  zinc-blende  — depending  on  which  .side  we  are  of  the  order-di.sorder  transilioii  line—  solid  soliitioii  plia.se 
a-s  tlierniodynamic  eqnilibriniii  state.  Tlie  corresponding  phase  would  llins  be  frozen  in  at  low  temperature.  The 
argnnietil  behind  this  is  that  at  temperatures  lo*v  with  respect  to  the  melting  point  phase  .separalioii  would  be 
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Ijiiuleretl  liy  kiii»‘tic  harri<*r'  which  we  will  se-  are  larje  In  a  «iigh(  variation  of  this  picture,  one  imagines  the  orfler- 
cli.sorder  transition  line  continued  inside  the  miscibility  gap  as  a  metastable  transition.  In  that  case,  if  we  are  at  a 
point  vertically  below  a  thermodynamic  et^uilibrium  liiic-blende  phase,  but  on  the  left  of  the  metastable  traitsiiion 
line,  the  diamond  phase  would  have  a  lower  free  energy  than  the  zinc-blende  phase,  even  though  both  have  higher 
free  energies  than  the  phase  separated  state.  The  argument  that  such  a  transition  could  occur  is  that  the  change 
from  zinc-b|emle  to  diamond  ordering  requires  only  local  rearangement  of  the  atoms.  This  is  less  hiiidererl  by  kinetic 
Itarriers  than  is  pha.se  .separation.  Of  course,  the  speed  at  which  the  temperature  of  the  system  is  lowerwl  coiihl 
determine  which  of  the  the  two  outcomes  prevails. 

We  note  that  there  has  been  significant  controversy  regarding  this  picture  in  eiiiier  of  its  two  variants.  A  pro|'ose,| 
alternative  theory  is  that  the  order-dtsorder  transition  between  the  diamond  and  zinc-blende  phases,  found  experi¬ 
mentally  for  the  Cie-ba.setl  alloys,  is  purely  of  surface  growth  kinetic  origin  [21.22].  .Still  another  possiliilily  snggestetl 
In'  Osorio  et  al.  (2()]  but  remaining  to  be  e.xplored  is  that  a  surface  thermodynamic  etinilibriinii  phase  grts  frozen  in 
becau.se  of  absence  of  full  bulk  thermodynamic  et|uiiibratioii. 

Otie  way  to  test  the  validity  of  the  "bulk  thermodynamic'  approach  to  this  problem  is  to  calculate  the  relevant 
interaction  parameters  from  first-principles.  There  have  only  been  limited  attempts  to  calcnlate  the  interaction 
parameters  of  the  5  =  1  Hamiltonian  from  electronic  structure  calculations.  The  e.xLsling  calculations  [2'l-26j  leatl  ‘o 
parameters  whicli  would  give  a  very  high  miscibility  temperature  (significantly  higher  than  the  liqiiiihi.s  line).  Those 
restdfs  are  similar  to  our  present  findings  for  the  C’-BN  system.  Osorio  et  al.  (2t5)  foimil  this  to  be  true  even  when 
AZ  =  ±2  boiuls  were  e.xcluded  and  the  electrostatic  and  AZ  =  ±1  compen.sation  effects  were  taken  into  accoiml 
as.  in  effect,  has  been  done  in  the  present  work  by  our  a-ssumptions  about  short-range  order.  Now.  a  nii.scibilily  gap 
teiiiperatnre  above  the  melting  point,  as  depicted  in  Fig.  4(a).  makes  the  .scenario  ofqnenching  in  the  zinc-Weiide  form 
from  a  higli  temperature  solid  solution  impossible.  This  fact  lends  indirect  support  to  the  idea  that  surface  (kinetic 
or  p.sendo-ei|uilibrium  thermodynamic)  effects  need  to  be  invoked  to  e.xplain  the  zinc-blende  ordering.  Intleetl.  a 
miscibility  gap  temperature  above  the  melting  temperature  would  indicate  that  the  nii.xed  crystals  or  alloys  can  only 
be  formed  by  conden.sation  or  solidification  of  a  mixed  gas  or  liquid  phase  in  which,  of  course,  a  zinc-blende  ordering 
does  not  make  sen.se. 

The  “.solid  soliilions"  of  c-BN-diamond  obtained  in  Badzian's  work  [5]  were  obtained  by  direct  (i.e.  nou-catalytir ) 
solid  state  transformation  of  he.xagonal  C-BN  crystals  (obtained  by  CVD)  at  high  temperature  >  33W  K  ami  pressnr»-s 
of  —  14  GPa.  As  stated  in  the  paper,  temperatures  above  the  melting  temperature  may  have  occurerl  in  the  center  of 
the  sample.  Thus,  the  synthesis  of  mixed  cry.$tals  does  not  contradict  our  prediction  of  the  rnkscibility  gap  letn|>eraHire 
being  above  the  melting  point. 

In  summary,  by  analogy  to  the  more  extensively  studied  Ge-ba.sed  ternaries,  the  pha.se  diagram  of  Fig.  4(a) 
following  directly  from  onr  calculations,  appears  more  likely  than  that  of  Fig.  4(b).  Nevertheless,  we  note  lltai  none 
of  the  calculations  of  lieterovalent  ternaries  .so  far  has  included  the  type  of  long-range  effects  which  were  fonnd  to 
l>e  important  for  common-anion  p.seudobiiiaries.  We  can  thus  not  completely  e.xchide  the  situation  iiulicaletl  in  Fi^. 
4(l>).  The  experimental  evidence  is  at  present  insufficient  to  distinguish  between  the  two  pictures. 

Reinxiiig  .some  of  the  a-ssumptioiis  about  short-range  order  in  our  treatment  (e  g.  the  requirement  that  B  anil  N 
alw.nys  occur  as  nearest  iieigiibors)  woidd  introduce  configuralioiis  of  higher  energy  and  thus  presumaWy  lurnasr 
T\ii;.  Only  the  existence  of  Special  long-range  order  wonhl  .seem  to  allow  for  an  appreciable  lowering  ofTin;.  Even 
with  a  lowering  of  Tmc  by  6U  %.  it  does  not  lie  much  lower  than  the  liquidns.  The  ronchision  of  limited  niLvihilily 
at  rea.sonably  low  temperature'  would  tlius  remain  valid  even  in  this  ca.se. 

Althongli  the  C-BN  alloys  are  thus  “in  principle'  unstable  towards  plia.se  separation  al  room  temperature,  it  ran 
eri.<;ily  be  that  alloys  are  hindered  from  pha.<=e  separating  by  large  kinetic  barriers.  Indeed,  we  note  that  phta-v*- 
.■'••pjiral  ion  would  require  boiifl  breaking  and  sul^staiilial  diffusion.  Since  the  bond-'trengtlis  are  of  the  onler  of  4  e\'. 
kiiieiir  barriers  will  be  huge  compared  to  the  "driving  force"  towards  pha.se  separation.  The  latter  can  l>e  iileniilh-tl 
witli  the  energy  of  formation  which  is  an  order  of  magnitude  smaller  as  mentioned  above. 

For  all  practical  purposes,  these  mi.xed  crystals  thus  appear  like  true  solid  solutions  even  Ihoiigli  they  are  in  principle 
nietastalile  As  is  well  known,  this  .system  is  in  addition  metastable  towards  the  hexagonal  layered  plia.ses  at  amluiMit 
pressure. 


C.  Properties  of  disordered  alloys 

The  properties  of  the  di.sordered  alloys  .such  as  the  energy  of  formation  and  the  baiulgap  ran  l>e  e.xpanderl  in  a 
cluster  e.xpansion.  We  note  that  within  the  CA'M  one  might  calculate  the  probability  distribution  of  the  telrnhetlral 
clii.sters  from  the  energy  minimization  of  the  free  energy.  These  may  deviate  from  the  random  (hinoiuial)  ilistribniioii. 

However,  since  the  miscibility  gap  critical  temperature  is  rather  nnrertaiii.  as  was  di.vnssed  above,  and  llie  nii.xe.1 
componmis  are.  in  fact,  likely  to  have  been  formetl  in  the  liquid  state,  the  ii.se  of  a  random  disiribniion  for  stal«lir,al 
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H  li»*r»*  (li«*  i*ii<*rgy  of  formation  of  the  orderetl  rompound?  AJ/..lj)  vra.s  suKslitiited  from  Eq  (2).  This  staiistirni 
averaging  is  seen  to  lead  simply  to  a  renormalization  Hy  a  factor  of  3/4. 

For  t  he  handgaps  of  the  disordered  alloys,  we  obtain  in  a  similar  fa.shion. 
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where  ’  is  the  handgap  of  the  n-th  CW-structure  asisociated  with  the  corespondiiig  tetrahedral  chister.  As  experted. 
the  resulting  band-gap  bowing,  shown  in  Fig.  1  by  the  full  continuous  line,  while  still  bring  vey  strong,  is  somewhat 
less  pronounced  than  that  of  the  ordered  compounds.  We  note  that  the  tetrahedron  appro.\iniatioii  itseil  for  our 
cluster  e.'cpansion  is  probably  less  converged  for  the  ba.’.dgaps  than  for  the  total  energies. 


V.  CO.NCLl'SION 

The  main  finding  of  this  work  is  that  c-BN-diamoiid  mi.Ked  crystals  and  alloys  e.'chibit  an  anomalously  large  hand-gap 
bowing,  which  we  propose  to  be  related  to  the  staggered  band-offsets  at  the  corresponding  helerojunction. 

We  also  di3cu.s.sed  the  phase  stability  within  the  tetrahedron  appro.vimation  of  a  cluster  e.'cpansion.  Our  iikhW 
of  the  alloys  is  based  on  the  absence  of  B-B  and  N-N  bonds  (because  of  their  e.'cpected  high  energy)  and  nearest 
neighlior  clustering  of  B-C  and  C-N  bonds  to  ensure  local  charge  neutrality.  Within  this  restrictive  model  for 
the  short-range  order,  the  problem  can  be  treated  as  a  p.seudobiuary  alloy  problem.  The  energies  of  formation 
within  this  model  were  found  to  map  onto  the  ferromagnetic  Ising  model.  Their  values  lead  to  the  prediction  that 
the  miscibility  gap  temperature  would  be  well  above  the  liquidus  line.  Although  this  conclusion  is  in  agre<Mitetit 
with  more  complete  treatments  of  other  heterovalent  ternary  .semiconductor  alloy  systems,  some  uncertainty  remaiiws 
with  regards  to  this  conclu.sion.  Indeed,  long-range  correlations  have  not  been  taken  into  account  in  any  of  these 
treatment.'*.  These  correlations  have  been  found  to  lower  the  miscibility  gap  temperature  by  as  much  as  60  %  in  .soni'* 
runimon-atiion  semicomluctor  alloy  sy.sleni.s.  An  alternative  tentative  phase  diagram  was  con.Mrnrietl  inrhiiling  .*nrli 
a  retuinii.ilization.  The  e.X|>erimen(al  conditions  in  which  the  mi.xed  crystals  were  synthesize*!  conhl  be  coiusisteni  with 
either  of  the  tentative  phase  diagrams  of  Fig.  -1  although  our  an8ly.si.s  seems  to  favor  Fig.  4(a).  Further  work  will  be 
riMpdreil  to  sort  out  which  of  these  two  alternatives  applies. 

The  “in  principle"  mela<«tahilily  of  the  system  should,  however,  not  be  a  practical  impediment  to  their  n.se  ami 
furl  her  study  since  kinetic  barriers  towards  phase  separation  are  e.'cpected  to  be  huge,  even  at  f.'irly  elevale*l  tein|«*r- 
atnres.  We  hope  this  study  will  lead  to  further  e.xperimental  work  on  this  system  which  conhl  verify  onr  predirtinns 
al'ont  the  hand  gaps. 
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FIG.  1  Band  gap  as  a  function  of  composition  for 
zinr-blende  derived  ( B.N )jC'j( j _j.,  ordered  crystals:  LDA 
(. '<111,1  res),  self-energy  corrected  [13]  (diamonds)  and  random 
allot  '  ('inuolli  full  line). 


FKi.  .’.  Band-ofl’sel  at  diamond-C'/f'Bi'l  helerojunction. 


FIG.  3.  Free  energ.v  of  (BN  )jC'j|  i_,|  alloys  within  the  reg¬ 
ular  .solution  model. 


FK  !.  4.  Tentative  phase  diagrams  of  C-BN-diamond.  (a) 
Ferrumagnelit'  tsing"  niutlel  treated  within  regular  solnlion 


ap|iroxiiiia(ioii  renormalized  by  92  %  lb)  same  w  la)  buc  i»- 
diidinR  an  ad~hor  60  %  lowering  of  the  laiedbility  gap  d«e 
to  lonfz-range  correlaiione.  A  liypotheiiticai  diamond  (0)  to 
zinr-bleiuk-  |ZB)  order-disorder  iraasition  is  tentatively  indi- 
catetl  with  ?  and  dotted  line.  The  Uqiiklus  line  separating 
the  litinid  |L)  from  the  solid  portion  of  the  phase  diagram  b 
indiraietl  in  both. 
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